Pull counting statistics for the Kondo dot in the unitary limit 
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We calculate the charge transfer probability distribution function xW for the Kondo dot in the 
strong coupling limit within the framework of the Nozieres-Fermi-liquid theory of the Kondo effect. 
At zero temperature, the ratio of the moments Cn of the charge distribution to the backscattering 
current 7bs follows a universal law C„/2/bs ~ ( — !)"(! + 2")/6. The functional form of xW is 
consistent with tunnelling of electrons and, possibly, electron pairs. We then discuss the cross-over 
behaviour of x(-^) from weak to strong Coulomb repulsion in the underlying Anderson impurity 
model and relate this to the existing results. Finally, we extend our analysis to the case of finite 
temperatures. 
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PACS numbers: 72.10.Fk, 71.10.Pm, 73.63.-b 

The importance of the noise spectra in charge trans- 
port has been recognised long ago [j]. The higher mo- 
ments of the charge distribution, however, are only just 
becoming accessible experimentally |2|. The probability 
P{Q) of charge Q being transmitted through the system 
during the measuring time T and the respective moment 
generating function x(A) = J2q e^^'^P{Q) have been es- 
tablished by Levitov and Lesovik J3l for non-interacting 
electrons. The effects of electron-electron correlations 
on the counting statistics are currently a subject of in- 
tensive debate. In particular, a universal result that the 
linear response, zero-temperature statistics is binomial 
for any interacting system with one conducting channel 
has recently emerged Q. 

An important test system of such type is the so-called 
Kondo dot: a semiconductor quantum dot device in the 
Kondo regime. Extensive experimental [5| and theoret- 
ical [a, Q inquiries into conducting properties of these 
systems have recently taken place. In the last few weeks, 
two works appeared, Ref. jl] and Ref. ^, where the zero- 
temperature shot noise is calculated for such devices in 
the strong-coupling limit. 

In this Letter we widen the investigation and calculate 
the full charge distribution function (which encodes all 
the moments) in the strong-coupling Kondo regime. We 
also discuss the temperature corrections and make con- 
tact with the previous results for the Anderson impurity 
model |1|. 

A good place to start is indeed the Hamiltonian for a 
single level quantum dot (which is a version of the An- 
derson impurity model) 



H — Hq + if T + Hq 



(1) 



where 



^0 = ^(ep-AiK)<,,«^p,,,« + 5^(Ao + fT/i)4d, , (2) 



P,I7,K 



describes the conducting leads, i^l^a,K. being the creation 
operator for the electron with momentum p, spin cr ='\ , [, 



in the lead k ~ R/L, with dispersion relation Sp such that 
the Fermi level density of states i' = ^ S{^p) i^ finite (we 
set Ep — and e = fi, = fcs = 1 throughout), and the 
fermionic level dj. at the energy Aq in magnetic field h. 
The leads arc DC biased ^l ^ — ^^ = V/2 > 0. Further, 



^T = E (^i^-^' 



'^a,L 



■ lRii,Rda 



H.c 



) (3) 



represents the local electron tunnelling between the leads 
and the dot (with, in general, different amplitudes 7l and 
7fl) while Hq = Un-^n^ with ria = d\d^, stands for the 
Coulomb repulsion U . In equilibrium, this probably is 
the best studied model in the condensed matter theory. 
In particular, in the Kondo regime (Ao < 0, [/ positive 
and large) , the Schricffer-Wolf type transformation [l^ , 
tailored to the lead geometry p , can be applied to result 
into a Kondo type or s-d exchange model. 

For the Kondo model in the strong-coupling limit, 
when the spin on the dot is absorbed into the Fermi sea 
forming a singlet, Nozieres [Ij devised a Landau- Fermi- 
liquid description based on a 'molecular field' expansion 
of the phase shift of the s-wave electrons: 



5a{e) = 5Q + ae + (l>''{n„ - Ug) 



(4) 



where 5q — 7r/2, a and (jf" are phenomenological param- 
eters corresponding to the residual potential scattering 
and the residual interactions, respectively. These pro- 
cesses are generated by polarising the Kondo singlet and 
so are of the order ~ l/TV , where Tk is the Kondo tem- 
perature. The specific heat coefficient is proportional to 
a/{T:u) while the magnetic susceptibility is proportional 
to the sum a/( -Kv) + 20°/7r. Simple arguments were ad- 
vanced in Ref. [Ill to the effect that, because the Kondo 
singularity is tied up to the Fermi level, there exists a 
relation a = 2v(j)'^ between the two processes in Eq. Q. 
In particular, this explains why the Wilson ratio is equal 
to 2. Clearly, it is not difficult to write down a sec- 
ond quantised Hamiltonian describing the processes in 



Eq. Q). Before doing so let us remark that it is well 
known that there is no critical U in the Anderson model 
(P) and the Fermi liquid approach is actually valid for all 
U (see [i3 for a review). 

The strong-coupling Hamiltonian that describes the 
scattering and interaction processes encoded in Nozieres 
Eq. (gj is of the form H = Hq + H^c + i?int- The free 
Hamiltonian here is 



Hn 



p.a 



■y{cl.aCp,a + «p,<Tap,<T) + yQ 



(5) 



where c) is the creation operator for the s-wave electrons, 
a^ is the creation operator of the p-wave electrons, in- 
cluded in order to account for the transport |7| , and the 
operator 



<3= T^'Y^icl.aap^a 



r) 



stands for the (minus) charge transferred across the junc- 
tion. The scattering term is 



H,. 



E(^ 






' ^P' J^p.a^p'^a 



(6) 



while the interaction term reads 



gauged away by the canonical transformation 



c - 
a 



c\ = cos(A/4)c — isin(A/4)a , 
■ a\ = — 2sin(A/4)c + cos(A/4)a 



(8) 



We therefore reach an important conclusion that the 
charge measuring field enters this problem as a rotation 
of the strong-coupling basis of the s- and the p-states. 
While Hq is invariant under this substitution, it should 
be performed in both the scattering and the interaction 
Hamiltonians, iJgc [c] + i?int [c] ^ Hx = Hsc[cx] + Hint[cx], 
when calculating the statistics. It is easily checked that 
at the first order in A: Hx = ffsc+i?int + (A/4)/bs+0(A2), 
where /bs is the backscattering current operator 



^(4a, - aUMc^ , (9) 



Atti'Tk 



2tiv'^Tk 
alternatively available from the commutator /bi 

The charge probability distribution function is there- 
fore defined by 



X(A) 



AN\ 



T, 



cexp 



Hx{t)dt 



c 



(10) 



Hi-nt = 



ffl t t 



(7) 



where c„ = ^ c„ ^ and we have changed to the dimen- 
sionless amplitudes a and (f>, so that in the actual Kondo 
model a = (J3 = 1 (in the intermediate calculations it is 
convenient to treat a and as free parameters though). 
By the nature of the strong-coupling fixed point, the 
operators a and (p are irrelevant in the renormalisation 
group sense and therefore the perturbative expansion in 
a and (j) is expected to converge. 

The charge measuring field A(t) [A(t) = X9{t)6{T — t) 
on the forward path and A(t) = —X9{t)6{T — t) on the 
backward pass of the Keldysh contour C], couples, in the 
Lagrangian formulation, to the current |13J via a term in 
the action J dtX{t)Q{t) = — J dt\{t)Q{t), which can be 



where Tc is the time ordering operator on the Keldysh 
contour, and N = {VT)/n = {2e^/h)VT. Applying 
the standard linked cluster expansion (still valid on the 
Keldysh contour, of course) [l^, we see that the leading 
correction to the distribution function is given by a con- 
nected average 



Inx(A) = iNX 



The neglected terms a^ 



dtidt2{Tc{Hxitl)Hx{t2)})c + - 

(11) 

(jp, 0^, etc., are of the higher 
order in voltage (temperature) than the main correction 
because of the irrelevant nature of the perturbation. In 
order to make progress with Eq. (|ll|l . one only needs the 
Green's function of the A-rotated c-operator, which is 
easily seen to be the following matrix in Keldysh space 



J 



gx{p,u;) = nrSiep - lu) { [/{u; - V/2) + f{iu - V/2) - 1]% + [er'^'^fiuj - V/2) + e'^'^f{uj + V/2)]t+ 

-[(1 - fiuj - y/2))e'V2 + (1 _ j(^ + F/2))e-'V2]^_ | ^ 



(12) 



r 



where fi is the standard choice of Pauli matrices and /(w) 
is the Fermi distribution function. We did not write the 
principal part which does not contribute to local quanti- 
ties in the flat band model [la. 



The correction to the distribution function which is 
due to the scattering term © is shown in Fig.^a). The 
corresponding analytic expression is (a factor of 2 comes 
from summing over the spin index) 



(^alnxlA) 






2 p 



dcc;tj^[(e-*^ - 1)71l(1 - n^) + (e^^ - 1)^^(1 - "i)] 



(13) 



which, at zero temperature, contributes to Eq. Hll() a 
term 



5ah\x{X) 



12nT^ 



(e- 



1) 



(14) 



Regarding the correction to the charge distribution 
coming from the interaction term Q , any diagrams with 
a single insertion of the Green's function vanish (therefore 
there is also no a4> cross term) and the only remaining 
connected graph is shown in Fig. ^b). This is best cal- 
culated in real time. At zero temperature we have (no 
summation over the spin index): 



S^lnxW = 



dt 



2n^v^Tl jc 



dhdt2g{h,t2f9{t2,hf 



^'Tl 






cos^[X/2 + {Vt)/2)] 



e-'^ - 1 






(15) 



(e 



-2iA 



1) 



Combining the results we find that the zero- 
temperature charge distribution function is of the form: 



Inx(A) 



F^T 



iNX 

2(j)\e-'^'^ ~ I)] +0{V^ 



UttTJ. 



[(a2 + 0^)(e-*'-l) 



(16) 



In particular, the results of [2| for the average backscat- 
tering current and noise power are correctly reproduced 
(there the notation is a = a/(27r),6 = ^/ti")- Therefore 
the Fano factor is equal to 



C2 
2/bs 



2(a2+502) 



(17) 
0) and 



which is 1/2 for non-interacting electrons (< 
5/6 at the Kondo fixed point {a ~ (p). 

But we now know much more - all the mo- 
ments can be computed from Eq. p6|) via C„ = 
(-2)"(d"/dA")lnx(A)U=o, leading to 



-^ = (-1)^ 
2/bs ^ ' 



(1 + 2"+i) 



2(a2 + 502) 



(18) 



n = 2,3,.... Hence at a = = 1 we have a hierar- 
chy of universal (Fano factor inspired) ratios C„/2/bs = 
(-l)"(l-h2")/6. 



(a) 





FIG. 1: Schematic representation of the scattering (a), and 
interaction contributions (b). 



In some systems, like in the case of the fractional quan- 
tum Hall effect, there are physical grounds on which to 
expect fractional charge quasiparticles to appear in trans- 
port [13 • In the present problem, which after all answers 
the Fermi liquid description, there are no such circum- 
stances. In fact, the functional form of the distribution 
function allows a very simple interpretation. The term 
proportional to V^ {a^ + (jP'){e^'^^ ~1) must be interpreted 
as tunnelling of conventional electrons with charge e (this 
is indeed obvious from the fact that at = there is no 
interaction and the problem is trivial). It is tempting to 
interpret the term proportional to V^(j)'^{e^'^'^^ — \), which 
appears as a result of electron correlations, as a coher- 
ent tunnelling of electron pairs with charge 2e. [Certain 
caution is required here, as one can show that, at higher 
orders, a term of the form V^ a^ {e~'^'^^ — 1) exists, which 
is an artefact of the expansion around the perfect trans- 
mission.] 

Let us now return to the original model Eq. J^l- There 
is an extensive vintage literature on the Fermi liquid 
properties of the equilibrium Anderson model [I3 . The 
perturbative expansion in the Coulomb repulsion [/ can 
often be re-summed to all orders, using Ward identities 
and similar techniques, so that the observable quantities 
(specific heat, conductance, magnetisation) are expressed 
in terms of the even at odd susceptibilities Xe and Xo (cor- 
relations of n| with n| and rij , respectively) . These are 
universal functions of V and are known exactly from the 
Bethe ansatz [13, [13 • Specifically, in the weak coupling 
;t, = 1-)- (3 - 7r74)C/V(7r2r2) + ... and Xo = -(7/(7rr), 
where F = 2(77^/7)^ (we assume, for simplicity, a symmet- 
ric junction, 7^ = 7^ = 7 and particle-hole symmetry 
Aq = — t/). On the other hand, in the strong coupling 
we have: Xe = {Va)l{'KTK) and Xo = {V^)I{i^Tk), as 
here a = = 1 and TV is the Kondo temperature up 
to a prefactor [13, [la [IS ■ The programme of extending 
a Fermi liquid approach to non-equilibrium properties of 
the Anderson model has not been comprehensively car- 



ricd out yet. 

There is a Fermi- liquid proof, due to Oguri |2Cll |. 
that the leading non-equilibrium correction to the zero- 
temperature current is of the form 



/bs = 



yi 



127r2r2 



(Xe + ^Xl) 



(19) 



which is valid for all U and interpolates between the 
weak-coupling and the strong-coupling regimes of the 
Anderson model. We see that the above result for /bs is 
simply the strong-coupling limit of Oguri's formula. 

As to the noise and higher moments no analogous 
Fermi-liquid results exist, to the best of our knowledge. 
However, in our recent paper [J|, we did put forward the 
formula for the charge distribution function: 



Inx(A) 



N liX + 



3r2 



Xe 



xl 



ie 



-iX 



(20) 



Xo / - 



2iX 



)]} 



which was a guess based on the properties of the per- 
turbative expansion in U. In jj| we have shown that 
Eq. I|2()(l is correct at the order U^ [it is also consistent 
with Eq. H19|) for all [/]. The calculations in this Letter 
prove that Eq. H20|l is also correct in the strong-coupling 
limit. 



Before we close, we wish to discuss the effects a finite 
temperature T. The energy integral in Eq. (|13(l is elemen- 
tary. The time integration in Eq. (|15|l can also be done 
with standard methods (the real time Green's function 
at T > is related to that at T = by a well known 
conformal mapping 1/t — * ttT/ sinh(7rri) , see |2l|). The 
result for the generating function is: 



Inx(A) = iNX + 

202 



."' +f^j T/^ + 4(7rT)2][e^/2T(e-'^ - 1) + e-^/2^(e*^ - 1)] 



VT 



247rr2. |_sinh(y/2T) 



sinh(y/r) 



[1/2 + (^r)2][e^/^(e-2a _ ^) ^ g-y/T(g-2.A _ ^^j 



(21) 



Using this formula we recover the finite-T version of Eq. (|19|l . see [23|. The higher moments generated from Eq. 121|) . 
starting with n = 2, are all new results (to our knowledge). The most important of these is the expression for the 
noise, which we here report: 



V ( ( V 

Co = TT { coth — 



12 



-V^ + — — {irTf 



coth(;)f [y2 + (.T)2]} 



(22) 



r 



To summarise, we have calculated the full counting 
statistics near the unitary strong-coupling limit of the 
Kondo dot. We wish to stress that Eq. H2(J|) . even though 
correct for all the moments in both weak-coupling and 
strong-coupling limits, has not been formally proven yet. 



The theoretical challenge is to develop Fermi liquid the- 
ory for the counting statistics. 
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